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ABSTRACT 


Boiling Heat Transfer with Low Steam Quality 
Thomas Frank Marucci 
Submitted to the Department of Naval Architecture and Marine Engineer- 
ing on May 21, 1965, in partial fullfillment of the requirements for 
tne Master of Science Degree in Mechanical Engineering and the Profos- 
sional Degree, Naval Engineer. 


Heat transfer for water and steam mixtures in a horizontal stain- 
less steel tube under non-boiling, surface-boiling, and fully-developed 
boiling conditions was investigated. 


Local heat transfer coefficients are presented for a limited range 
of flow parameters. Inside diameter of 0.2125 inch and a range of 1h.) 
to 78.5 L/D are considered. Exit pressure varied from 60 to about 72 
psia for mass velocities from 5.9 x 102, to 2.38 x 106 lbs/hr ft. The 
range of steam quality covered was O to lh%. The majority of this data 
was age at inlet temperatures of about 2909F. Heat fluxes ranged from 
2 x 10? to 1.8 x 106 Btu/nr ft2 and were limited by burnout conditions, 


Heat transfer results for non-boiling and subcooled nucleate boil- 
ing are in agreement with previous investigations. Bulk boiling data 
is presented at the two lower flow rates. This data shows that as quali- 
ty is increased, the local heat transfer coefficient also increases. This 
is indicated on the boiling curve as a shift toward lower wall superheat 
as quality increases. Similar results are indicated to a lesser degree 
with subcooled liquids at low values of subcooling. This appears to be 
due to the relatively large volume fraction of vapor created in the nu- 
cleate boiling process. 


Experimental heat transfer coefficients obtained are considered to 
be in fair agreement with predicted values. 


Thesis Supervisor: Arthur E. Bergles 
Title: Assistant Professor of Mechanical Engineering 
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Surface area 

test section cross section area 
specific heat 

diameter 

voltage across the test section 
shunt voltage 

mass velocity 

enthalpy 

heat transfer coefficient 
current 

thermal conductivity 

heated length 

pressure 

heat transfer rate 

heat flux 

heat flux at the inception of boiling 
critical heat flux 
forced=convection heat flux 


electrical resistance 
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Nu 
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Re 


temperature 

inside wall temperature 
outside wall temperature 
velocity j 

mass rate of flow 


position along the test section measured from the 
inlet to the heated section 


steam quality 

temperature drop through the tube wall 
wall - minus - saturation temperature 
saturation - minus - bulk temperature 
dynamic viscosity 

surface tension 

density 

resistivity 


latent heat of vaporization 


Dimensionless Groups 
Nusselt number = hD/k 
Prandtl number ^ C, AA [x 


Reynolds number = CD/AM 





BO 


exit 


sat 


Subscripts 


bulk or mixed-mean conditions 
burnout condition 

inlet bulk condition 

test section exit bulk condition 
liquid condition 

saturation condition 

tube 

vapor condition 


condition at heat transfer surface 
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INTRODUCTION 

Ine current development of high performance heat-transfer 
equipment nas stimulated a great amount of research in boiling 
heat transfer. Past studies of boiling heat transfer have been 
a factor in the development of nuclear reactors and rocket motors. 
Vaporization of the coolant is frequently desired, such as in the 
boiling-water reactor and bootstrap rocket engine. 

There has been a vast amount of work in the field on forced 
convection-boiling heat transfer. Most of this has been concen- 
trated on the prediction of burnout under given conditions of 
fluid, surface, and flow parameters. These studies have produced 
numerous correlations; the best of which can predict critical heat 
flux to about plus or minus 30$ of experimental data. 

Within the past ten years there has been an increase in the 
number of investigations designed to study heat transfer rates en- 
countered in boiling heat transfer. The results of these studies 
have been presented in two basic forms. First, in the form of the 
boiling curve where heat flux is plotted versus wall superheat 
1-1), This has Usually been utilized where the bulk fluid is 


subcooled, and indicates that heat flux can be represented as: 
n 
αλ οσα -ᾱ ο (1) 


Other studies considered bulk boiling and present results that show 
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the local heat transfer coefficient (h) increasing as steam quality 
increases beyond a certain range. This range has been generally 
agreed upon as beginning at about 15% quality. As quality varies 
below this value h is constant. 

Generally, when liquid flows through a heated tube, subcooled 
nucleate boiling can take place before the bulk of the fluid reaches 
its saturation point, if the impressed heat flux is high enough. The 
equilibrium thermodynamic quality would be zero in this subcooled 
region of flow, but the void fraction could be substantial, depending 
of course on the fluid properties and heat transfer rate. As flow 
proceeds down the tube the bulk of the fluid becomes saturated, net 
vapor generation commences. In this portion of the flow, the vapor 
generation can take place by additional nucleate boiling and/or by 
evaporation at the liquid-vapor interface. If the tube is long enough, 
it is possible to reach the point of 100% vapor. 

This study was conducted to extend the fully-developed boiling 
curve in the region of bulk boiling and to determine what effect of 
variation in h with quality has on the boiling curve. It was also 
pointed out by Professor Bergles that non equilibrium voids created 
in subcooled region should have an effect on h. 

The equipment used limits the study to low pressure, low flow 
rates and moderate heat fluxes. Fluid conditions included low sub- 
cooling and low steam quality. This study is particularly applicable 


in the development of evaporators for saline water conversion. 
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PROCEDURE 

Operation and Jata Taking 

Ihe experimental apparatus (Figure 1) used in this study was 
One already in operation at the M.I.T. Heat Transfer Laboratory, 
Only the test sections had to be fabricated to suit this study. 
Generally, the test sections had to be instrumented to provide 
means of determining local pressure and the temperature at the 
heat transfer surface. To obtain these measurements, pressure 
taps were located along the tube and thermocouples were placed 
on the outside wall of the test section. To insure that the out- 
Side wall was adiabatic, two guard heaters were used and main- 
tained at a temperature equal to the tube wall temperature. A 
brief discussion of the experimental apparatus and a detailed 
discription of the construction of the test section is presented 
in Appendix A, A diagram of the heat transfer test section is 
shown in Figure 2. A preliminary investigation of critical heat 
flux was conducted to be able to predict the critical heat flux 
for the range of parameters to be studied. This was done pri- 
marily to avoid burnout of the elaborately instrumented heat 
transfer test sections. The test sections used to determine 
critical heat flux were not instrumented to obtain wall temp- 
eratures and were constructed with only one pressure tap located 
at the test section exit. 


Prior to installation, the test section was given a final 
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internal cleaning with acetone. After installation, the entire test 
section was wrapped with insulation to minimize external heat losses. 
The loop was then filled with distilled water and all lines vented to 
remove air. The degassing operation was then initiated by heating the 
water in the degassing tank. With the water boiling, the loop water 
was heated and slowly sprayed into the top of tne tank. Results from 
Reference (1) indicated that 20 minutes of such degassing reduced the 
dissolvec-air content to approximately 1 c.c. /liter. Two phase flow 
can me luence system stability and recent studies have indicated that 
system instability can cause premature burnout. Air trapped in the 
system would provide a certain amount of compressibility. Reference (2) 
discusses system instabilities in more detail. After degassing, the 
heat exchanger cooling water was turned on to bring the system to nor- 
mal temperature. The loop was then ready for operation, circulating 
distilled, deionized, and degassed water. 

Critical heat flux data was obtained oy increasing power gradu- 
ally while maintaining constant flow rate, exit pressure, and inlet 
temperature. Critical heat flux was evaluated at the last recorded 
conditions orior to physical burnout of the test section. To vary 
exit quality conditions, inlet temperature was varied in subsequent 
runs. 
| The general procedure for taking heat transfer data was similar 
to that used for critical heat flux data. The major difference was that 


at each power level, while other conditions were maintained, the guard- 
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heater power was adjusted to the balance point, i.e., where the average 
guard-neater temperature (two-thermocouples) was the same as the aver- 
age outer tube wall temperature (three thermocouples). These tempera~ 
tures could not be kept exactly equal, but were maintained to within 
plus or minus 5°F of each other. At equilibrium, wall temperatures, 
fluid temperatures, pressures, pressure differences and test section 
voltage and current were recorded. The heat flux was then increased 

to the next level and the process repeated. 

After the highest heat flux was completed, the power was reduced 
and the inlet temperature changed. Readings were then taken for a set 
of heat fluxes at a new inlet temperature. 

The maximum attainable heat flux for any given run was determined 
by burnout except in one instance due to voltage limitations of the 
generators. The burnout heat flux was estimated from data presented 
in Reference (3) and from data taken in this study. It was desired 
to approach but not reach turnout. Temperature and pressure levels 
were largely determined by the system. Beginning at the lowest inlet 
temperature attainable, the inlet temperature was increased to about 
20°F below the saturation temperature for the corresponding pressure. 
The exit pressure was maintained at 60 psia except at final runs where 
the minimum exit pressure obtainable for the desired flow rate exceeded 
60 psia. This was attributed to the fact that choked flow conditions 
were present due to the large volume of vapor content in the test section. 

Initial data taken with T.S.-1 and 2 was at relatively low values 


of heat flux and at low inlet temperatures. This was to obtain non- 
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boiling heat transfer data which is well correlated and hence serves 

as & check on instrumentation. By increasing the heat flux slowly it 
was possible to obtain data on the inception of nucleate boiling and to 
be able to trace the boiling curve in the transition region between non- 
boiling and fully-developed subcooled nucleate boiling. 

At the higher heat fluxes with intermediate inlet temperatures the 
fully-developed subcooled nucleate boiling curve could be well estab- 
lished. 

The final data was then taken with T.S.-2 at inlet temperatures 
near saturation. The inlet temperature was maintained below saturation 
at all times since the equipment did not allow determination of inlet 
conditions when net vapor existed at the inlet. This procedure pro- 
vided data for bulk boiling over the quality range from 0 to 21%. To 
extend the data to higher steam qualities it was necessary to fabri- 
cate T.S.-3. 3T.S.-3 was constructed similar to T.8,-1 and 2 except 
that it was twice as long. The location of pressure taps and wall 
thermocouples for each test section is given in Appendix A. Due to 
inconsistent results as discussed later, T.5,-3 was modified by moving 
the thermocouples at the upstream locations to positions corresponding 
to L/D's at the downstream position of T.S.-2. One series of runs was 
conducted with this test section at G=1.19 x 100 and er 290°F 


Accidental burnout and time limitations precluded going further with 


this study. 
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Data Reduction 

All data reduction was done by hand calculation with the aid of a 
desk calculator, This included computation of local flow and heat trans- 
fer parameters, including flow rate, heat flux, wall superheat, bulk temp- 
erature or enthalpy as well as the dimensionless groupings Re, Nu, Pr, 
Nu/PrOoh, The mass velocity (G) was determined from the fluid flow rate 
and tube dimensions. The temperatures were converted from thermocouple 
emfs with NBS standard calibration. The heat input was obtained from 
the tube voltage and current (EI) and checked by the resistance and cur- 
rent, a third value was obtained from the mass flow, specific heat, and 
temperature rise. Values obtained by each were in excellent agreement 
and differed at most by 10%, but were usually within 5% of each other. 

The latter method was of no value when the exit temperature reached sat- 
uration. The average heat flux was calculated using the several values 
of heat input. 

Heat transfer results were based on meaSurements of the outer wall 
temperature. 4 correction for temperature drop through the tube was made 
to first obtain the inner wall temperature. A modified form of the Krieth 
and Summerfield solution was used to obtain the necessary correction. 

This solution assumes an infinite cylindrical resistor, with heat trans- 
fer at the inner surface and an adiabatic outer surface. The final form 
of the equations as well as thermal-property information is given in 
Reference (1). The final form of the curve(„T,, versus I at various values 
of Tow iS shown in (Figure 3). It should be noted that at high values 


of current, the wall temperature drop is extremely large. The data taken 





ο 

in the last run (Table IV Series N Appendix C) points out the effects 
of this on the results. The wall temperature drop is so large that the 
inside wall temperature is actually reduced as the power level is in- 
creased. Results then indicate a substantial reduction in the degree 
of wall superneat. This leads to the conclusion that the equation for 
wall temperature drop may not produce accurate results at high values 
of current. High currents were reached with T.S.-2 as well, (Table IV 
Series F) but with a shorter test section the outside wall temperatures 
were also higher, hence, the effect was not as pronounced, At each 
thermocouple position, calculations were made to obtain the local heat- 
transfer results. The fluid bulk temperatures are based on a linear 
temperature gradient where data is for subcooled bulk fluid. Where net 
vapor existed,the bulk temperature is the saturation temperature for the 
particular local pressure. The local pressure was obtained graphically 
by drawing a smooth curve through the observed pressures at three loca- 
tions. All fluid properties are based on these local bulk temperatures» 
With the saturation temperature (T at), the degree of subcooling CATS 
Er. E ») and the wall superheat (AT..4 = τ. = Tat) were available. 

A sample set of calculations is shown in Appendix B. A tabulated 


summary of the data obtained is given in Appendix C, 
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RESULTS 
Critical Heat Flux 
Burnout data were required to set an upper limit on the heat flux 
that could safely be applied to the heat transfer section. This data 
is presented graphically in Figure h, for G = 1.19 x 10° 1b/hr ft? 
(V ~5ft/sec). Also plotted are data taken by Wessel (3)* for the same 
tube size with varying fluid velocity and a slightly larger tube with 
varying L/D(V ~ 5ft/sec). The data shows an increase in (yA)... as 
exit quality increases. 
Heat Transfer Results 
The non-boiling results of this study are presented in Figure 5, 
along with the line indicating McAdams (lh) correlating equation. Data 
in Figure were obtained for G = 2.38, 1.19 and 0.59 x 106 lb/hr γος 
(Υ 510, 5, and 2.5 ft/sec). The span of Reynolds number also results 
from the bulk temperature variation. The data is consistently above 


that predicted by McAdams equation and can be represented by: 
Nu = 0.017 Re? 85 PrOel (2) 


Subcooled nucleate boiling data is plotted in Figure 6. Figure 6 
contains data at the same flow rates, P ~ 60 psia, L/D = 37.4, and a 
range of subcooling from about O to 175%F, Three other curves are 
plotted to aid in describing the results. 

First, the curve of (α/λ}ῃρ is the curve of heat flux versus wall 
superheat aue) obtained from non-boiling prediction. This curve 


is based onG = 1,19 x 10° lb/hr ft?, απ. 109-183°F, where 


* Numbers denote reference listed in Bibliography 





¿52 


/A E m P _m _ ; l 

(α, dec A(T n) πας, T d D T) and h is determined from 

properties at bulk temperature corresponding to the data using Eq (2). 
Second, the curve of (q/A), represents the locus of points at which 


Surface boiling begins and where the heat flux departs from the non- 


boiling prediction. Tnis is discussed in Reference (5) and is given by: 


0.0234 


1.156 2.30/p 


τ... 


W sat 


(a/A), = 15.60 p (3) 
The curve plotted here is for p = 60 psia. 

The third curve represents the mean line througn data of Wessel 
for fully-developed subcooled nucleate boiling and was determined as 
indicated in Figure (7). 

The data covers a range of flow conditions in the transition region 
from the inception of boiling to a point on the fully-develoned boil- 
ing curve at which bulk boiling begins. As fully-developed boiling is 
reached, the data approaches an asymptote which agrees with ‘vessel's 
data. Figure (6) illustrates that in the transition region heat trans- 
fer rates are deoendent on voth flow rate and subcooling, and as the 
fully-developed boiling region is approached the dependence on flow 
rate and subcooling disappears. 

Incipient boiling does appear to begin at the point predicted by 
the intersection of the incipient boiling curve and the non-boiling 
prediction. 

Bulk boiling data is presented in Figures 8 and 9. Figure 8 shows 


data for L/D = 37.4, at the two lower flow rates for bulk boiling 
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(i.e., at points where the bulk fluid was at the saturation tempera- 
ture). The data represents the range of steam quality from O to 21%, 
Ine fully-developed curve is again plotted for comparison. The data 
of Figure 8 to appear to follow the fully-developed boiling curve for 
the range of quality indicated. 

Figure 9 contains data obtained with T.S. -2, 3, and 3 modified. 
The range of steam quality is extended to almost 10%, and data is 
plotted for varying L/D's. The data for a particular L/D can still 
be represented by a form of Eq. (1). But, as quality increases (rep- 
resented by increasing L/D), the fully-developed boiling curves are 
shifted toward the lower wall superheat. The effect of quality seems 
to override the pressure effect, i.e., the lower pressure results in 
higher wall superheat. 

Similar results appear in Figure 10 for data with subcooled 
boiling conditions. As L/D increases, at first, wall superheat in- 
creases, and then decreases. This appears to be independent of the 
degree of subcooling within the range covered in Figure 10 ( ον» 
30 150%), 

In an attempt to determine the effect of steam quality on heat 
transfer coefficient, Figures ll through ll show the experimentally de- 
termined heat transfer coefficient versus steam quality. Figure ll con- 
tains data obtained with T.S. -2 for G = 0.59 and 1.19 x 10° lb,/hr rt’. 
The results indicate that heat transfer coefficients are dependent on 


heat flux, and flow rate alone in the range of steam quality. At 
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quantities above 10-15$ the data appears to merge for each flow rate 
and the dependence upon heat flux disappears. The effect of flow rate 
is quite pronounced as would be expected. At lower flow rates the effect 
of quality on the local transfer coefficient becomes apparent at lower 
values of quality. 

Figure ll shows the results of T.S. -2 and 3 at G = 1.19% 10°, and 
Figure 12 shows similar results for G = 0,59 x 10°, Data obtained from 
these test sections appears to be inconsistent. Results of T.S. -2 in- 
dicate an increase inh as flow proceeds downstream. If the data ob- 
tained from T,5. -3 is then compared with tnis, it would appear that h 
is relatively constant between 37.4 L/D and 57.9 L/D. Beyond that point 
h reaches a maximum and then decreases. Because of the inconsistency in 
the data, T.S. -3 was modified. The thermocouples located upstream on 
T.S. -3 (L/D = 57.9 and 64.1) were moved to positions which corresponded 
to the downstream thermocouple positions of T.S. -2 (L/D = 30.9 and 37.4). 

Data taken with this test section is shown in Figure 13. Also plot» 
ted for comparision is data taken with T.S. -2 and 3 at G = 1.19 x 10° 
lbs/hr ft’. Data for each test section was taken at approximately equal 
values of heat flux and inlet temperature, in so far as data was avail- 
able from the other test sections. Results of these last tests compare 
favorably with that obtained from test sections 2 and 3. The absolute 
values of h at corresponding L/D's are about the same. The only differ- 
ence is that data from the modified test section shows an increase inh 


at the last thermocouple position instead of the decrease noted with 


179. -3 
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In order to compare the bulk boiling data with previous work in 
this area, Figure li shows the dimensionless ratio h/hy versus the 
Martinelli parameter πο where the Martinelli parameter is defined 
Dy: 


ΝΕ (x/1-x)0-? T Zu 


AR (hy 


h is the experimental heat transfer coefficient and h is the liquid heat 
trandfer coefficient as given by the McAdams correlation. Also plotted 
here is Chen's F function and his prediction of h at parameters approx- 
imating the mean values for the data. The curve lies below the data at 
the lower L/D's of T.S. -2 , and 3 modified, but above the data for high 
L/D's of T.S. -3 and 3 modified. Chen's correlation is presented in a 
later section. 

Figure 15 shows similar data at G = .59 x 10° along with Chen's F 
function and his predition of h at parameters approximating the mean 
values for the data. Here Chen's F factor seems to divide the data for 
low L/D's and high L/D's. At the low flow rate Chen's prediction of h 


is higher than the experimental values. 
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DISCUSSION OF RESULTS 


Critical Heat Flux 

In the course of obtaining this data several interesting phenomena 
occurred, Although the primary intent of this study was to obtain heat 
transfer data, it is considered necessary to mention these items briefly. 

First, the fact that (a/A) increased as much as 50% as the exit 
quality increased from 0 to about 25%. Previous studies in this area 
with subcooled fluid conditions at the exit have indicated that (yA), 
decreases as the amount of exit subcooling decreases. The most recent 
studies show that as zero subcooling is approached (9/4). appears to 
level off and in some instances increases . This is indicated by the 
data of Wessel shown in Figure kh. One of the best known studies of 
burnout with exit quality is that of Lowdermilk (6). His results in- 
dicate that E is at a maximum when exit quality is approximately 
20$ depending on other flow parameters as well. The data did not cover 
the low range of exit quality. Hence, previous work indicated that, at 
à constant pressure and flow rate,(o/A)... decreased as the exit subcool- 
ing decreased. Critical heat flux reaches a relative minimum, increases 
in the low range of exit quality to a relative maximum, and then falls 
off sharply. Bergles (1) also shows an increase in critical heat flux 
with an increase in quality up to 12%. Hence, the data presented here 
is in agreement with other studies. This critical heat flux is not as 
high as that attained with subcooled exit conditions. For more infor- 
mation in this area the reader is referred to recent studies by Loosmore 


and Skinner (7) on critical heat flux with subcooled boíling and by 
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Lopina (8) on critical heat flux with exit quality conditions. 

Second, the actual location of burnout of the test section was 
not always at the exit. In fact, at exit quality near zero, and with 
exit quality about 25%, burnout did occur at the exit. Between these 
limits, the location of burnout was considerably upstream of the exit. 
In one run with exit quality approximately 13.6% burnout occurred at 
18.7 L/D, while the overall tube dimensions in all runs was πα ΠΠ} 
Calculations indicate that the bulk fluid condition at the burnout 
location was slightly subcooled. This indicates that in the low range 
of exit qualities, exit conditions may not be the correct parameter to 
consider in determining critical heat flux. Lopina (8) and Waters, etal, 
(9) nave also observed upstream burnout and discuss this in more detail, 
Waters! study shows that upstream burnout is reproduceable and not neces- 
Sarily caused by instability. 

Third, it was observed during operation that as tbe exit quality 
increascd, the pressure drop from the test section exit to the exit 
plenum chamber became quite large. Hence, when measuring critical 
heat flux with net vapor generation, exit pressure must be determined 
at the test section exit and not in an exit plenum chamber. Results 
of this study indicated that there was essentially no pressure dif- 
ference between these two locations at zero quality, but at 17.8% 
quality the pressure difference was approximately 12.5 psia at G = 
19x 10° lbs/hr ft?. To explain this it is necessary to give a 
brief discription of the method used to control the pressure. 


The exit pressure level was controlled by means of a plug valve 
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downstream of the exit. At low values of heat flux, this val 
about midway between the fully opened and closed positions in order 
to maintain the desired pressure. As the heat flux was increased, 
and boiling conditions were reached, the exit pressure began to in- 
crease. This valve was then opened to reduce the pressure to the 
desired level. When bulk boiling conditions were reached, it was 
observed that pressure control was no longer possible, since the 
valve was wide open. As indicated previously, when the longer test 
section was used, pressure control was lost prior to obtaining the 
maximum heat flux for that series of runs. This then illustrates the 
choking effect that causes the pressure level to increase in the test 
section. 

As the steam-water mixture leaves the test section, it enters a 
plenum chamber made up of 3/4 inch brass pipe. Expansion of the steam 
can take place immediately. The pressure drop encountered from the 
test section exit and the plenum chamber is the result of this rapid 
expansion. 

Heat Transfer Results 

Non-boiling data shown in Figure 5 plots consistently above the 
correlation of McAdams. The divergence appears to be primarily due to 
the radial property variation which is not properly accounted for in 
the correlation. This temperature difference effect is clearly demon- 
strated in Reference (10). In any case, Eq. (1) is not intended to be 
a generally valid correlation, but is merely indicated to represent the 
data. 


Figure 6 shows the portion of the boiling curve from the inception 





203 
of boiling to the fully-developed boiling curve for three flow rates 
and two ranges of subcoolings. As indicated previously, inception of 
boiling did occur at the predicted value of wall Superheat and heat 
ουκ. Beyond this point, the slope of the curve increases until fully- 
developed-boiling is reached. Using the data of Wessel to establish 
the fully-developed boiling curve (Figure 7), the data does approach 
this line as an asymptote. Hence, the subcooled data appears to be con- 
sistent with the work of previous studies. 

At the onset of this study, Professor Bergles pointed out, that at 
high heat fluxes with low subcooling, tne data appeared to be shifted 
from the fully-developed curve, toward lower wall superheats. It has 
already been mentioned,that fully-developed subcooled boiling data 
(Figure 10) show a change in wall superheat at different positions a- 
long the length of the tube. At the inlet, subcooled liquid velocity 
can be specified as V. As the fluid moves down the heated sections, 
surface boiling begins and non-equilibrium voids are formed. Therefore, 
the volume occupied by the liquid is reduced. In order that continuity 
be maintained, the liquid velocity must increase. Let this new velocity 
be Vo, As boiling becomes more vigorous along the tube, more volume is 
Occupied by vapor and the liquid velocity continues to increase, This 
illustrates how the liquid volocity may increase along the length of the 
tube. With this in mind,and noting from Figure 6 that by increasing 
fluid velocity, fully-developed boiling is reached at higher values of 
heat flux. 

If we consider three fluid velocities which do not differ by much, 


and determine the heat flux at which the lowest velocity reaches the 
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fully-developed boiling curve, it is clear that at this same heat 
flux, the wall superheat will be reduced at the two higker velocities. 
This example shows how the wall superheat can decrease at higher L/D's 
with subcooled liquids. This also indicates that the local heat trans- 
fer coefficient must increase in this range. The heat transfer co- 


efficient is defined as: 


pa) een) (5) 


As long as the bulk fluid is subgooled, the bulk temperature is in- 
creasing along the direction of flow, If (Ty - Toat) decreases as 
indicated previously and Ίσας decreases with pressure, then Ty must 
also be decreasing. As a result, (1, - T.) must decrease and the heat 
transfer coefficient increases due to the presence of void fractions. 
This explains the trend toward lower wall superheat in data at high 
L/D's as shown in Figure 10. 

Bulk boiling data (Figure 9) also shows a reduction in wall super- 
heat with increasing L/D. Bennett, etal. (11) presents a discussion on 
the effect of net vapor generation, similar to the discussion above on 
the effects of void fraction on fluid velocity. This explains the lower 
wall superheats obtained at the exit of the test section. In the region 
of bulk boiling wnere the bulk temperature is the saturation temperature, 
à decrease in wall superheat means an increase in h for a specified heat 
flux. Therefore, wherever net vapor exists, there is an increase in the 
local heat transfer coefficient. 

Now, consider a section along the tube where the bulk fluid is slight- 


ly subcooled. As flow proceeds down tne tube the bulk fluid approaches 
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saturation conditions. The rate at which bubbles nucleate at the 
surface decreases since the fluid velocity is greater. At the same 
time, the bubbles, which have just departed the surface, do not 
collapse and the volume fraction of vapor remains nearly constant. 
Then, the fluid velocity will also remain constant and the degree 

of wall superheat predicted by the boiling curve will not change. 
Therefore, in this region near zero quality, the heat transfer co- 
efficient is relatively constant as well. Further along the heated 
section, net vapor generation continues by means of surface boiling 
and/or evaporation at the liouid vapor interface. At this point, the 
fluid velocity begins to increase again and continues to increase 
until the test section exit is reached. This appears to explain the 
data of Figure 9, as well as that of Figure 11. 

Figure 12 with data of T.S. -2 and 3 at G = 1.19 x 10° los/hr fi? 
shows the discrepancy in the data obtained from these two test sections. 
Data obtained with T.S. -3 shows that, as x increases, h rcaches a max- 
imum and then decreases. This could not be explained adequately. Con- 
sidering data for an approximate equal heat flux for the two test 
sections, there also appears to be a discontinuity in the data between 
L/D = 39.4 and 57.9. The value of h increases as L/D increases to ο 
but apparently there is no further increase in h to L/D = 57.9. Beyond 
this point, h reaches its maximum value and decreases, Results for G = 
$.9 x 105 lbs/hr ft? shown in Figure 12 are similar. 

It was expected that data obtained with T.S. -2 modified would re- 


' solve this discrepancy. Figure 13 shows data for one flow rate and one 
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value of heat flux for T,S. -2, 3, and 3 modified. The results of this 
show that the absolute magnitudes of h remained about the same. Της only 
difference which did appear consistently was that at high L/Dts. nume 
creased rather than decreased. The heat transfer coefficient then show- 
ed and increase over the first half of the T.S., then remained a rela- 
tively constant value, and finally increased again near the exit. One 
possible explanation is; at the lower L/D's flow conditions had just 
reached the point at which the void fraction is nearly constant. Over 
a short range beyond this point, h is constant with x. If this were 
true, the results as indicated in Figure 11, are misleading and quality 
does not significantly affect h until it reaches a value of about 20%. 

Another method of evaluating the data would be to draw the mean 
line through the data from both test section. It is uncertain which 
data, if any, is actually incorrect. If this were done in Figure 12, 
it would indicate the effect of quality becomes apparent at about 15$, 
but its overall effect would not be as great as indicated in Figure 11, 

Figures 15 and 16 show the same data on another set of coordinates, 
h/h, versus l/ X, e Dengler and Addoms (12) were tne first to use these 
coordinates in evaluating this type of data. Here Dr is defined as the 
conventional liquid heat transfer coefficient determined by McAdams! 
correlation. They proposed a correlation for two phase heat trensfer 
coefficients in terms of Lae, (the Martinelli parameter) and an F 
factor to account for the degree of nucleate boiling present. As the 
fluid progresses into the annular flow regime, with continuiously in- 
creasing steam quality, the increased velocity of the two phase nixtttm 
induced by the vaporization process suppresses the nucleate boiling proc- 
ess. Beyond this point the heat transfer coefficient becomes governed 


by the forced convection process. 





ips 
With this in mind Chen (13) proposed the following correlation 


for the two phase heat transfer coefficient: 


h = hnic * hpac (6) 


where mic is the microconvective portion due to the boiling mech- 


anism and hence dependent on liquid and bubble parameters: 
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and id is the macroconvective given by: 


er - nx = (8) 


T 


Again hy is given by McAdams' equation. F is a function of the 
Martinelli parameter and is defined as the ratio of the two phase 
Reynolds number to the liquid Reynolds number. S is called the 
Suppression factor and approaches unity at zero flow rete, and zero 
at infinite flow rate. Gouse (1) suggests that this suppression 
factor should go to zero much sooner and thus acknowledge the possi- 
bility of total suppression of nucleation at finite Reynolds number. 
Both S and F are presented graphically in Reference (13). 


Chen restricts the applicability of this correlation to systems 


with annuler or annular mist flow in vertical geometries. This 
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correlation successfully predicts the data of numerous independent 
studies to within plus or minus 15%, and seems to be one of the best 
available. Gouse (1h) also states that there is no reason why it 
should not apply to annular flow in horizontal tubes. 

Data obtained with thermocouples at the top and bottom of the 
tube wall, at the same location along the length, produced no indica- 
tion of stratification. Thermocouples indicated that the temperature 
difference obtained in all cases was less than the accuracy of the 
measurement. Also, the greater temperature fluctuated between the top 
and bottom. With this in mind Chen's correlation is considered a 
Suitable means of checking the data. 

As noted earlier Chen's prediction seems to plot between the data 
for high and low L/D's indicating that the difference can be considered 
as scatter in the data. It must be noted, that scatter of such a con- 


sistent nature is unusual, 
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CONCLUSIONS 


Critical Hest Flux: 


The critical heat flux increased as much as 50% as exit quality 


was increased from 0 to 25%, 


meat Transfer 


l. 


The local heat transfer coefficient, obtained at subcooled 
boiling conditions, does appear to increase due to the pre- 
sence of voids formed during the boiling process. 

Ás bulk fluid approaches saturation conditions the local heat 
transfer coefficient is constant over a low range of steam 
quality. The extent of this range increases as the flow rate 
increases. 

For the flow rates considered in this investigation, the local 
heat transfer coefficient shows a significant increase as qual- 
ity is increased beyond 15%. 

Chen's correlation appears to be fairly accurate in the low 
range of steam quality, but it's accuracy improves as quality 


increases and nucleate boiling is suppressed. 
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RECOMMENDATIONS 


X 


Critical Heat Flux 


n. 


When bulk fluid conditions at the exit contain vapor, the 


exit conditions must be determined by the pressure at the 

test section exit rather than the pressure in an exit plenum 
chamber. 

Since upstream burnout is vossible with exit quality conditions, 
and stable flow, it may be of interest to determine if some para- 
meter, other than exit quality might give a more meaningful in- 


dication of the results. 


Heat Transfer 


l. 


Additional experiments should be conducted with test sections 
similar to T,S. -3 and instrumentation covering a greater por- 
tion of the test section. This is suggested to see if there 
is actually a region along the test section.where the local 
heat transfer coefficient is constant. 

Additional study of heat transfer coefficient in the low qual- 
ity range is still required, in order to develope a more accu- 
rate prediction in the range where heat transfer rates are de- 
pendent on the combined effects of boiling and convective heat 


transfer. 
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APPENDIX A 


taj 


APERTMENTAL APPARATUS 
Sa PASS D 


Flow Loop and Pover Supply 

Ihe flow loop used in the experimental program was one already in 
operation at ihe M.I.T. Heat Transfer Laboratory. This facility is 
shown schematically in Figure 1. It was a closed loop system in which 
all components are of corrosion resistant materials, The system con- 
tained the main circulating pump, an accumulator, the test section line 
with parallel flow meters, oreheater, test section and its instrumenta- 
tion, the bypass line, and a heat exchanger utilizing city water. Aux- 
iliary equipment included a fill pump, supply tank, degassing tank, and 
a continuous deminer&lizer. The test section power was supplied by an 
a.C. motor - d.c. generator set. A discription of the major components 
follows. More detailed information is given in Reference (1). 

The main circulating pumo was a two-stage, turbine-type pump driven 
by a 3hp induction motor. It provided a head of 250 osi at 3.6 gpm. A 
bladder type accumulator, nressurized with nitrogen, located at the pump 
exit served to damp out pressure fluctuations. Tne ball valve in the by- 
pass line was used to control the by-pass flow rate and pressure, wnich 
in turn determined the pressure in tne test-section line. The test sec- 
tion line contained a Fischer-Porter flow meter and a set of preheaters 
in series. There are four 5kw preheaters plus an additional 5kw emersion 
unit fabricated in a 3 inch copper tube and provided with 0 to 100% con- 
trol by a Powerstat auto-transformer. This arrangement allowed varying 
preheat from O to 25kw, one being variable and the other allowing incre- 


ments of 5kw each. 
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Ihe test section power was supplied by motor-cenerators. The 
generators were driven by 440 volt, 3 phase synenronous motors. Two 
36kw d.c. generators, each nominally rated at 12 volts and 3000 amps 
were connected in series. The generators vere provided with water 
cooled shunts in parallel with the test section which allowed them to 


ve run open circuited at the test section for starting, or after burn- 


out. 
Instrumentation 


Instrumentation was available for reading pressure levels, differ- 
ential pressurc, temperature, test section flow rate and test section 
voltage and current. Pressure levels were read on bourdon-tube gages 
as shown in Figure 1. At tne test section inlet and exit were 200 and 
100 psi test gages with an accuracy of plus or minus 1/h% of full scale. 

Test section flow rate was determined from the larger flowrater meter 
which had interchangaole tubes and floats. The ranges covered are 20- 
160, 35-00, and 140-1360 lb/ir. áll units were calibrated as installed 
Moe system. 

The power supplied to the test section was obtainea from the test 
section voltage and current. The voltage was read directly on a weston 
multiple-range d.c. voltmeter with a specified accuracy of plus or minus 
1/2%. The current was determined from a standard shunt with a calibra- 
tion of 60.17 amp/mv. 

All temperatures were measured with copoer-constantan thermocouples 


made from 30 gage duplex wire. The fluid bulk temperature at the inlet 








p 
and exit of the test section wss measured by inserting thermocouples 
directly into the fluid using conax fittings with lava sealants. 
Methods of measuring wall temperature are deccrived later. The out- 
put voltages of the thermocouples and the shunt were displayed on a 
continuous recorder. The recorder was a Brown, single channel in- 
strument having ranges of 0-5, 5-11, 10-16, 15-21, and 20-26 mv. 

The manometer system used to measure differential pressure con- 
sists of two Meriam 60 inch U-tube manometers, manifolds and valves, 
and connecting lines of rubber nose. The system is designed to read 
a maximum of 10 pressure differences with either one or bod reference 
points on either manometer. One manometer was filled with mercury 
which gave a maximum range of ‘approximately 25 psi. The other mano- 
meter contained an oil with a specific gravity of 2.00 which allowed 
a maximum oressure dirference of approximately 2 psi. A typical test 
section with three pressure taps is showm in Figure 2. The pressure 
drop from the reference "A" located at the inlet, to points 2, and 3 
could be rezd on either the mercury or oil manometers deperding on the 
magnitude of the pressure drop. The two pressure gages were connected 
to the inlet and exit to give the pressure level and € check on the 
pressure drop. Vents were included so that all lines could be purged 
2 air. 

Measurement oi Tube-"all zemperature 
The most difficult problem in any heat transfer experiment is the 


measurement of the temperature of the heat transfer surface. The basic 


method of measuring tube wall temperature is described in Reference (1). 
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The method used here was as follows: 

Heated snields constructed of standard aluminum tubing, with 
threadlike grooves cut into the outer wall, were wrapped with in- 
sulated Chromel A-26 wire. Spacers fabricated from aluminum served 
to hold the shield in place. 

Power was supplec in the Chromel heating coil by a 110 VAC supply 
which is stepped down to about 25 volts by a Stancor filament trans-~ 
former. The power to each shield is controlled by a Powerstat. This 
arrangement permits line control of shicld heating. 

A layer of Scotch Electrical tape Number 27 was wrapped around 
the tube at each thermocouple location. The measuring junction was 
then bound radially around the tube with another layer of tape.  Thermo- 
couples were also taped to the inside of the shield. Sufficient length 
of all thermocouples were left inside the shield to avoid conduction 
errors. The space between the shield and tube was filled with asbestos 
fiber. Precautions were taken to insure the shield was electrically 
insulated from the test section. 

Power to the shields is adjusted until the shield thermocouples 
and the tube wall thermocouples within the shield section indicate 
approximately ecual EMF. 

Test en. 

The test section was constructed from 30h stainless steel tubing 
with an inside diameter of 0.2125 inch. Brass bushings were cut from 
3/4 inch round stock. Holes were drilled for the test section. The 


bushings were then sweated to the tube with silver solder. Holes 0.099 
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inch were then drilled in the bushings for pressure taps. These 
holes did not penetrate the test section. A number 75 drill (0.021 
inch) was used to drill the pressure tap hole. Beryllium copper pres- 
sure tap tubes (2mm ID) were then set in the bushings and soft soldered 
in place. The inside of this tube was filled with milk of magnesia to 
prevent the tap being plugged. With this completed, the milk of magnesia 
could then be removed. 

Án additional pressure tap was placed near the center of the test 
section. Here a snort stainless steel tube was first filled with milk 
of magnesia and then silver brazed to the test section. Excess braze 
was then filed away to give a very small fillet at the test section. 

The milk of magnesia was removed and the pressure tap hole drilled (0.021 
inch) through the test section. Any burrs inside the tube were removed 
with fine emery cloth and pieces of steel wool pushed through the tube, 

Thermocouples were then attached to the test-section wall as des- 
cribed previously. Thermocouples were located at four positions along 
the length. At two locations the thermocouples were placed perpendicu- 
lar to the tube axis, and the lead wires wrapped around the tube and 
taped in place. At the other two locations, thermocouples were again 
olaced perpendicular to the tube axis, but at each of these positions 
two thermocouples were used, one on the top of the tube and the other 
at the bottom. This set up was intended to detect stratification. 

The guard heaters were then assembled onto the test section. After 
the first guard heater was in place a section of 0.099 inch OD tubing 
was placed over the center pressure tap and soft soldered to the stain- 


less steel tube to extend this pressure tap. A diagram of this test 
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T.S. -l and 2 were identical in construction. Accidental burnout 
early in this investigation is the reason for this. The location of 


pressure taps and wall thermocouples is indicated below. 


fee, 1 and 2 


Calming length 6 in. Distance to thermocouples 
Heated length 10 in. ἴον. τ HELD 
Center pressure tap at 27.3 L/D 2 23.1 

3 30.9 

L 39.4 
io. 3 
Calming length 6 in. 1 57.9 L/D 
Heated length 20 in. 2 oho 
Center pressure tap at 45.5 L/D 3 το η 

L το 


feo. > modified 


Calming length 6 in. 1 30.9 L/D 
Heated length 20 in. 2 39.4 
Center pressure tap at 45.5 L/D 3 ie 


L 78.5 











SAMPLE CALCU aon 


l. Calculation of heat input to the test section: 
(a) from the current and voltage measurements 


E os 3.113 EI ο θα ο ο... 


= 205.36 E E, Btu/hr 


S 


E and B were measured values, 


(>) from current and resistence 


qy * 3113 1? 3 - 1.236 xi0* z^ m 


ΠΔ μα, 
NIE: co and o e. ION 
&. - (374) (2? - o0^)Ahh = 2.115 x 107 σεξ 
R = 0394 x 10^ Po 
J was determined from Reference (1) as a function of the 


tube temperature. ine tube temperature is given by: 


T. mp (A 


D OW 


AT was determined from #igure 2 and is a function of To 
W 
I, Des was measured by means or thermocouples converted with 


N.B.S. tabies. 


(c) fron mass flow rate and bulk temperature change 


νυ So (T exit - 4) 


+ and T, as well as all other temperatures was obtained 


T : 
eX = 


as indicated above for ie 


A sample set of data and calculations are shown on the following pages. 





Data Cor Run A-8 


τ Bs Pi A ?i-2 | ic) P3 ms a 

V Ey psig An-oil in-ci psig mv F 
n c8 5.62 06.099. 029 6.9 HUSO Ου, ie 
q e 7 w μα. Te ne 
T 3 ο. Tou-lb m. low-2 
mv ay mv my F mv Sm 


Ü de 
Peon 11,50 11.37 δι 11.29 — 6. 11.40 162 


~ 


ate A a, 
ae 


- 


E T poet i 7 Ñ 
C D) QUESO. Tow-la tow=hb 
mv nv mv y ny 07 NV 07 


ES 11.53 311.66 170.7 11.51 1:56.7 dE S 


go, = 12.65 (18.62) (205.36) = 48122 Btu/hr 


( 
dy * 1.236 x 10? (18.62)? 1.119 x 107? « LTTSO Btu/ar 


I = 60.17 (18.62) = 1120 amps & Tg, = 465 OF 


+ 


τονε; 


T. = h65 - 57.5 = 107.5 °F and ja, = 2.84 


U 


2. eat transfer surface area A= TT(ID)L * 0.0529 ft^ 


ll. The change in enthaloy through the test section is 


DI g/w = 135 Jtu/1b, 


“Δ 1g, Tg and Tp denote guard heater terineretures. 


A 
4a 


resrectively. 


(1075) ohm ft 


Average nect flux a/a = 1:0093/ .0529 = 9.159 x 10? Btu/hr rue 


The change in bulk tempcrature from inlet to any position along the 


Tre a and b indicate temreratures at the top and bottom of the tube 





b (Texit ~ 71) x'/L 


The pressure at any point along the tube w.s obtained graphically 


ed 


as indicated earlier. 

Entering the stecm tables with the local pressure, the saturation 
temperature, tne saturated liouid enthaloy and the latent heat of 
vacorization were ootainecd. 


& 


all fluid properties at cach thermocougle position were obtained 


6 


us 


au the local bulk temperature. These include Co» k, P,4, and Pre 


At each position the following values were obtained as indicated 


acove. 
L/b x'/L D ANTE ^ P T 
: Y wW x Sub 
x! 
p.h .35 115  3h5.0 60.2 292.9 
Ὁ 58.7 


23.7 .575 162.0 115 347.9 59.9 292.6 


Eco — . 75 hyo.] X15 355.7 O 


37.4 .906 115 > 319.3 zo lg 


st a and b were averaged to give the outside wall 





L/D AT, D AT EE (1 
Ei le Dp 52.1 191.6 2/927 
2007 Τό 192 DOM NOS 215.0 
Bose 99 5o 6304 ον. 900.0 
3704 119.9 165.9 57.2. 126.2 163.4 


10. The local value of h was not calculated for subcoolec boiling data, 


but was determined for non-boiling and bulk boiling data as shown; 


A m ΕΝ 
di (ο/ 4) / Ce T) 


llo Nu and Re were determined for non-boiling data as well as the ratio 


Nu/PreeE | 


Nu = hD/k = (h/k) (0.0202) 
Re - GD/44 = 2.404 x 1 /4 @ w= 36h 1b, / hr 
Pr 9 OS M /k 


Pr is tabuleted as a function of temperature in Reference (15) with 


the other properties lor water. 


12. Bulk boiling data included some additional calculations as indicated: 


(a) 


(Ὁ) 


local enthaloy M nts = A 


9 


"here Ah! , 3 (hn! ,, - nop x!/L 


exit 


NOTE: When this value exceeded the value of saturated licuid 


enthalpy, net vapor was present. The local quality was 
calculated using the local enthalpy, saturated liquid 
enthalpy and latent heat of vaporization. 


Evaluation of the Martinelli parameter required values of AA 


and JO for steam. My was obtained from Reference Πώς 
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Py vas obtained from the steam tables. The Martinelli 
parameter was then calculated for each data point at which 
net vapor was present. The Martinelli parameter is defined 


as: 





MXS * CX? RNA AA 


To obtain the ratio h/hi, the liquid heat transfer coefficient 
(hr) was calculated at each of these data points with McAdams 
equation 

hy = 0.023 (x/D) ποθ prone 
This equation was used because it provided & means of compar- 


ing the data with Chen's work in this field. His value of hy... 


is equal to hy,as determined above, multiplied by his F factor. 


A summary of the data is contained in Appendix C, 
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SUNMARY OF SORTA 


Burnout Data 


Non-boiling Datz 


Subcooled Nucleate Boiling Data 


Bulk Boiling Data 


Range of Variables 


Page 


LO 


li 


LS 


55 


69 





D= 0,2125 inch 


RUN Pexit 


psia 
1 30 
2 56 
52 
LO 
1,2 
L5 
0.5.2 Iz 


Ov VU E LW 


i Burnout occured at tre exit 


oe 


TABLE I 


BURNOUT DATA 


L = 10 inches T, BM esc) EE 


A 2 X -= h! 
A E 

ie 325 20 

el 65 

2135 211 

1.18 52 

12500 109 

1.638 160 

1.95 296.1, 


V = 5 ft/sec 


(προ = D coc! 
Btu/lb 


K 


17 
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